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Graphene on a substrate has been shown to exhibit a transition, depending on the substrate material, from 
a zero-gap semiconductor state to a semimetallic state. The ground-state energy of the electron (hole) gas 
has been calculated within the random-phase approximation. 
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1. INTRODUCTION 

A monolayer of carbon atoms that form a perfect 
hexagonal lattice (graphene) has such a band struc- 
ture that the energy gap is equal to zero in three K 
and three K' points of the Brillouin zone [1]. The 
Schrodinger equation in the k • p approximation has 
the form [2] 

H F(r) = sF(r), (1) 

^=K' £ °*y (2) 

\ ucr ■ k I 

where the quantity u — §7a « 9.84 x 10 7 cm/s is 
similar to the Kane matrix element for the rate of in- 
terband transitions in the Dirac model [3], 7 ~ 3 eV 
is the band parameter numerically equal to the over- 
lap integral of atomic orbitals that participate in the 
chemical bonds of carbon atoms in graphene, ao = 
1.44 A is the interatomic distance in the graphene lat- 
tice, cr = (or, 02) are the Pauli matrices, k = (k x , k y ), 
k = {k x , ky) , k X y = idx,y: and H = 1. 

The unitary transformation U\ = ( p j 5 , ) , where / 
is the 2x2 unit matrix, reduces the Hamiltonian H 
to the form 

H' = UiH TJl = H* °X (3) 
\ —ucr ■ k I 

The corresponding equation with the transformed wave- 
function ip(r) = UiF(r) is equivalent to a pair of Weyl 
equations 1 . In quantum electrodynamics (QED), the 

lr The first application of the Weyl equation to the descrip- 
tion of the charge carriers in a zero-gap semiconductor was 
probably proposed in [6] in view of the appearance of the ax- 
ial Adler— Bell— Jackiw anomaly in parallel electric and magnetic 
fields. 



Weyl equation describes neutrino, a massless spin-1/2 
particle. Since the particles were assumed to be spin- 
less when deriving Eq. (1), Novoselov [1] introduced 
the notion of pseudospin. As correctly mentioned in 
[5], a complete model should include an 8 x 8 ma- 
trix Hamiltonian due to the twofold valley degeneracy, 
twofold pseudospin degeneracy, and twofold spin de- 
generacy. The 8x8 Hamiltonian may be reduced to 
the 4x4 Hamiltonian given by Eq. (3) if the Fermi 
momentum is determined in the spin-unpolarized state 
of the particles. In this case, the degeneracy multiplic- 
ity is v = v e h = 2 (z^e./i is the valley degeneracy of the 
conduction band or the valence band). 

To obtain the Dirac equation, let us perform an- 
other unitary transformation U2 = ^ ( \ lj ) 

HZ = U 2 H'M={ °r U(T n k )^u a -k, (4) 
\ucr • k I 

where a = ( ° ) are ^ ne Dirac matrices. Thus, the 
charge carriers in graphene are described in the frame- 
work of the zero-gap (A = 0) Dirac model. The use of 
the Dirac equation as a 4 x 4 matrix equation in the 
twodimensional system is justified because the 4x4 and 
2x2 matrix representations may be equally used in the 
case of two spatial dimensions [6] . This allows us to ap- 
ply the diagram technique of QED to the case of the 
two-dimensional system of Dirac fermions (graphene). 

In this work, it is shown that graphene on a sub- 
strate may undergo a transition from a zero-gap semi- 
conductor state to a semimetallic state, depending on 

the substrate material. The transition occurs at a* ~ 

2 

1, where a* = — — is the analog of the fine-structure 
constant, whose numerical value depends on the rel- 
ative permittivity e\ % °f the media surrounding the 
graphene and n e f f = £l jp 2 , similar to a thin film [7] . 



2. GROUND-STATE ENERGY 
(GENERAL CONSIDERATION) 

To determine whether the zero-gap semiconductor 
phase of graphene is stable with respect to the transi- 
tion to another phase, one has to calculate the ground- 
state energy of the electron (hole) gas that appears 
in graphene when the electric field is applied. The 
ground-state energy per particle is the sum of the three 
terms 



Eexch ~t~ &o 



(5) 



Here, the average kinetic energy is given by the expres- 
sion Ekin = \wpp because, according to Eq. (1), the 
dispersion relation of the charge carriers is linear near 
the K and K' points of the Brillouin zone: e p — ±u|p| 
(+ and — correspond to electrons and holes, respec- 
tively), pf = \J^niujv is the Fermi momentum, n 2 D 
is the areal density of particles, and v is the above- 
mentioned degeneracy multiplicity. If the Fermi level 
ep lies above e = 0, the charge carriers in the system 
are only the conduction electrons with a number of val- 
leys v e = 2; if ep < 0, then the charge carriers are only 
holes with Vh = 2. The position of the Fermi level may 
be changed by applying the electric field [1] . Both cases 
are obviously equivalent in the Dirac model. Below, we 
will consider the case of electrons for definiteness. 

The exchange energy is given by the diagram (see 
Fig. 1) 



E, 



exch 



d 2 pde d 2 kdu 
2n 2D J (2^) 3 (2tt) 3 



Sp{r"( P) e;k,w) 



x G (p, e) ^G (k, to)} 4°J (p - k, e - w) , (6) 

where the photon propagator is D$) (p — k, e — u>) ~ 
V (p — k) 8^8^ (we neglect the photon poles whose 
contributions to the integral with respect to the fre- 
quencies e and ui are on the order of (u/c) 2 ~ 1CP 5 
each, i.e., small compared to the contribution of the 
Green's function poles), V (q) 



is the Coulomb 



K e//|q| 

law in the two-dimensional case. 

The correlation energy is given by the formula 2 [8] 



1 



d kduj 



2n 2D J (2tt) 



x 



dX 



-A^(k)n 44 (Ma;) 
1 - Ai/U(k)n 44 (k,iw) 



■ A^y(k)n 44 (k, iu) 



(7) 

The polarization operator nj^k, itu) in the lowest 
order in the interaction (see Fig. 2) is calculated from 
the zero-approximation Green's functions [9] 



G (0) (P,£) 



up 



(e p - e + iS-)(e p + e - iS + ) 



(8) 



2 Formula (7) is derived in the nonrelativistic case and is the 
sum of the ring diagrams of all orders (as the most divergent dia- 
grams). The situation in the relativistic case is similar; therefore, 
the same formula with the corresponding polarization operators 
is used here. 




G(k,w) 

Fig. 1. Exchange diagram. 
G(°)( P , £ ) 




G<°}(p - k, e - w) 

Fig. 2. Zero-approximation polarization operator. 

G( P , £ ) 




G(p-ke-w) 
Fig. 3. Total polarization operator. 

where p = "f^p^, p^ = (p, ie/u) (pseudo-Euclidean met- 



rics), 8± — Ssign(ep ± e p ), S 
frequency integral [10] 



n^faiw) = 16 



(k- P ) 



d 2 P e 



-0. Evaluating the 
(\P\-PF) 



2e c 



(i 



(k-p 



2 ' 



(9) 



The total polarization operator II 4 4 (k, iu) given by the 
diagram in Fig. 3 is a renormalizcd quantity with re- 
spect to II44 (k, uo) due to the Coulomb interaction. 

3. RENORMALIZATION PROCEDURE 

The relation for the constant u = §700 was ob- 
tained analytically by linearizing the dispersion rela- 
tion in the vicinity of the K and K' points of the Bril- 
louin zone. The dispersion relation was, in turn, found 
by solving the Schrodinger equation with the Hamil- 
tonian in the tight-binding approximation. However, 
the interaction between particles as the interaction in 
a many-body system is disregarded in this approxima- 
tion. Thus, the constant u must be renormalized tak- 
ing into account the Coulomb interaction. The renor- 
malized value of u, generally speaking, depends on the 
electron density. 
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The problem is strictly formulated as a problem 
of solving the system of integral equations for the ex- 
act Green's function G(p,e), exact photon propagator 
D fJiU (q, f2), and vertex function T(p,e;p' ,e') (as well 
as for the polarization operator and self-energy). The 
approximate solution of the equations is possible, but 
seems too lengthy. We use here the renormalization 
theory developed in QED. In fact, the exact Green's 
function, photon propagator, and vertex function are 
replaced by the respective quantities in the lowest or- 
der in the interaction multiplied by the renormalizing 
constants. 

We proceed from two assumptions: 

(i) the charge coincides with the physically observed 
charge (ignoring the environment) 



eo; 



(10) 



(ii) the effective mass of particles remains equal to 
zero; i.e., the Coulomb interaction does not open a gap 
until the transition point to the semimetallic state 



0. 



(11) 



The gap opening is energetically unfavorable [11], 
which is confirmed in experiments. 

Let us quote the known relations of the renormal- 
ization theory 

r^ 1 = Zi7", (12) 

G=Z 2 - 1 G (0) , (13) 

D VU , = Z^D$. (14) 

According to the renormalization theory, the charge 
satisfies the relation [12] 



e = Z 1 1 Z 2 Z^ 2 eo- 



(15) 



Taking into account Eq. (10) and the Ward identity 
Z\ = Z 2 , we obtain Z 3 = 1, which was assumed in Eq. 
(6). The Green's function satisfies the relation 

G = G (0) + G (0) £G (16) 

with the formal solution 

G -i = G (o)-i_ S; (17) 

where G^ -*™ 1 = —up. Therefore, it should be expected 
that the inclusion of the interaction results in the renor- 
malization of u, the only parameter in the dispersion 
relation, taking into account Eq. (11). Looking for the 
self-energy in the form 3 

£(p,e) = Aup, (18) 

we obtain 

G-\p,e) = -u*p, (19) 

where u* = (1 + A)u is the renormalized u value. The 
self-energy is given by the expression (see Fig. 4) 

S(p,e) = i f d ^T^p-q,e-n;p,s)G(p-q,e-n) 

xyiVfon). (20) 

3 Below, we show that S is independent of the frequency e 
and S(p) = Aufp, which does not, however, influence Eq. (19), 
because po at p4 = ipo should be replaced by e/u*. 




G(p — q,e — 0) 
Fig. 4. Self-energy. 

Taking into account Eqs. (12) and (13) and the Ward 
identity, we come to 

s( P) e) = i / ^y i G(°)(p-q, e -n)7*'i>iS?(q,n). 

(21) 

After simple calculations (see below), we obtain 



A = a*I 



PF 



(22) 



I(x) 



tt/2 



111 



1 — x 2 sin tp I dtp, (23) 



i.e., the renormalized dispersion relation has the form 
e* p = ±u\p\(l + a*l(^y (24) 

Self-Energy Calculation 

Evaluating the pole integral with respect to f2 in 
Eq. (21), we obtain the expressions for the imaginary 
and real part, which will be considered separately: 



d 2 q uy 1 ^ - qi) 



im£(p) = - 
1 , 



) 2 

d 2 q wy 2 (p; 



2e p - q 



Q2). 



V(q)0(|p-q|-pf), 
(25) 



(2tt) 2 
d 2 q 



If 



V(q)6(\p-q\-p F ) 
2 V(q)0(\p-q\- PF ). (26) 



2 ' J (2?r) 

Let p be directed along q x , then pi = p and p 2 = 
0. Integral (25) is the difference between two integrals 
(with pi and q\ in their numerators), both diverging 
at q ^> p. Let us expand the integrand of the integral 
with qi in the scries to the ~ p/q 2 term, which gives 

a * 1 , 1c 
— U7 pin — , 
4 pf 

where the upper momentum cutoff q c ~ 2^/3-^/310 is 
introduced, which is on the order of the half-distance 
between the neighboring K and K' points of the Bril- 
louin zone (it is the momentum at which the linear 
dispersion relation breaks down). 

To evaluate the integral with p\, let us introduce 
the dimensionless variable z = g J p a ™ B f . By taking the 
integral with respect to z, the integral with respect to 
the angle is reduced to the form 



47T 



1*7 p / dtp In 
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where 



4. GROUND-STATE ENERGY 



JPf ~P 2 sm 2 tp _q c -pcosip 
p|sin.</?| p|smy| 

Taking into account that z c 3> 1, we may write 



U7V / dtp In (^z c + + z^j 



a 

47T 



_ 01 1 1 1c , ol 1 f y\sm<p\ 

~ -^r u 1 P ln H T" M 7 P / d ¥ ln — o — • 

2 p F An J 2 



Combining all these integrals together, we obtain the 
real part in the form 



ReS(p) = -U7 1 pln — - 

4 p F 



t/2 



-U7 p 



" r " hi I 2~ + 2~ 1 



sin 2 <p 



(27) 

The first term in Eq. (27) must be disregarded due to 
the following reasons: 

(i) q c 3> PFi and In ^ may be arbitrarily large at 
an arbitrarily small p f ; the renormalization coefficient 
for u becomes negative and, thereby, the further calcu- 
lation of the ground-state energy E gs is senseless; 

(ii) the corresponding contribution to the coefficient 
A is independent of the particle momentum; one may 
expect that the particles with p — > are insensitive to 
the renormalization of u and the particles with p —> Pf 
are the most sensitive to this renormalization. Thus, 
the condition A(0) — must be fulfilled. 

The integration with respect to the angle in the first 
integral in Eq. (26) yields zero; the non-zero contribu- 
tion is 



ZmE(p) = y 7°w<? c 



Up F , 



(28) 



where E(x) is the complete elliptic integral of the sec- 
ond kind. The first term is removed by the requirement 
ImZ Reg (0) = 0, i.e., 



IrriEitegip) 



a 

7T 











\Pf)_ 



l°up F , 



which represents the momentum-dependent shift of the 
frequency e in the renormalized Green's function 

G _1 (p, e) = -(1 + ^)«7P ~ * + ImZ Reg (p)} , 

However, the shift is a slowly varying function of the 
momentum and, therefore, may be replaced by its av- 
erage value, i.e., by the constant by which the inte- 
gration with respect to the frequency may be shifted 
when calculating the diagrams. Thus, the shift may be 
disregarded and we arrive at the result given by Eqs. 
(22) and (23). 



The average kinetic energy is now equal to 
2 

Ekin = ^Up F + 



it/2 



2a * f 2 , f , ( 1 1 A , • 2 \ 7 

— upf I x ax j In I - + - y 1 — x z sin <p J dip 

(29) 



0.0342a* 



up F , 



Hence, the contribution of the renormalization is < 
0.1; the respective contribution to the exchange and 
correlation energies is expected to be of the same order 
of magnitude (or even smaller). Therefore, to simplify 
further calculations, it reasonable to retain the linear 
form of the dispersion relation by averaging Eq. (22) 
in |p| and replacing u* with u* = (1 + A)u, where 

1 tt/2 

A = — / dx [ In ( — + — \J \ — x 2 sin 2 ipj dtp 



7T 











-0.0269a* 



(30) 



The renormalization of Eq. (6) for the exchange en 

ergy results in the equality E exc h = Z 2 ~ 1 E < fJ ch , when 
,(0) 



E exch is the exchange energy calculated with the non- 
renormalized Green's function, photon propagator, and 
vertex function: 



(0) 



where 

1 1 2ir 

J = J dx J dy J dx 

000 



a* J 
'~2tt 



■UPF, 



(31) 



(1 + cosx)xy 



y/x 2 + y 2 



2xy cos x ^ 



with — 0.915965 . . . being the Catalan's constant. 

The correlation energy may be calculated in the sec- 
ond order of the perturbation theory. The correspond- 
ing diagrams are renormalized by the factor Z^ 2 ', thus, 
E C orr — Z^ 2 Ecorr- 41 The further calculation with the 
use of the asymptotic expressions for the polarization 
operator in the lowest order in the interaction at low 
and high transferred momenta leads to the expression 
[11] 



e!2L 



1287T 



3tt 



25 
27 



27is J 



up F . 



(32) 



The final result for the ground-state energy reads 

8(0 + *) 



E n 



2tt - - 0.0342a* 
3 



3V2tt 



4 The number of renormalized vertices is half the number of 
renormalized Green's functions and the other half of the ver- 
tices remain nonrenormalized. Thus, the nth order diagram is 
renormalized by a factor of ■ 
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a* v /3tt 25 1 \ 

X 1 - 0.0269a* ~ 64^2^ \~8~ ~ 27 ~ 27^ J 

x 7 a " J J"^) 1/2 . (33) 

(1 - 0.0269a*) 2 J V v J y ' 

5. TRANSITION TO A SEMIMETALLIC STATE 
It is seen from Eq. (33) that the coefficient of 

1/2 r- 

u (usd/v) changes its sign at aj « 1.0204 for v = 2.° 
When this coefficient is negative, the creation of elec- 
tron- hole pairs in the system of the two-dimensional 
massless Dirac fermions becomes favorable taking into 
account that the absolute value of the ground-state en- 
ergy is a monotonically increasing function of the par- 
ticle density. This behavior is a manifestation of the 
system instability with respect to the Coulomb inter- 
action. In this phase transition occurs. 

As mentioned above, the gap opening is energeti- 
cally unfavorable; therefore, this is a transition from 
a zero-gap semiconductor state to a semimetallic one. 
Remarkably, the transition occurs depending on the a* 
value, which in turn depends on the substrate permit- 
tivity. Thus, the substrate material determines whether 
graphene has semiconductor or semimetallic proper- 
ties. It is also seen that the spin-unpolarized state 
with v — 2 is energetically more favorable at a posi- 
tive value of the ground-state energy, whereas the spin- 
polarized state with v = 1 is more favorable at the 
opposite sign. 6 As a result, the additional transition 
to the spin-polarized state occurs simultaneously with 
the transition to the semimetallic state. 

The closeness of the parameter to unity indicates 
that there is a close analogy between the case under 
consideration and the instability of the Coulomb field 
of the charge Z — 137 in QED with respect to the spon- 
taneous creation of electron-positron pairs (the fine- 
structure constant effectively approaches unity). It is 
possible that the exact value of a*, at which the de- 
scribed transition occurs, is actually equal to unity and 
the insignificant deviation of ap from unity obtained in 
this work is due to the inaccuracy of the approximation 
used. 

6. POSSIBLE EXPERIMENTAL OBSERVATIONS 

According to the experiments, the a* value depends 
on the substrate material. For example, K e ff — 5 and 
a* w 0.44 for the Si0 2 substrate and K e jf — 3 and 
a* w 0.73 [13] for the SiC substrate, but the condition 
a* ~ 1, which corresponds to K e ff — 2, is required for 
the transition. 

Formally, n e ff = 1 in a vacuum and graphene must 
be a semimetal. The overlap of the valence band and 

5 For v = 1, we obtain 1.0214 and o^i > a o> wmcn is 

important if the transition is approached from the a* < 1 side. 
In this case, the transition from the spin-unpolarized phase to 
the spin-polarized one occurs due to the change in the sign of 
E gs . 

6 A spin factor of 2 is already taken into account in the formula 
for pp. Thus, in our notation, v should be replaced by z/ e h/2 in 
the spin-polarized case without changing the formulas. 



conduction band in the semimetallic state may be es- 
timated as [11] 

6E~(b-~^up F , (34) 
where b = a*ap, which is only meaningful at a* > 

1 /2 

«q. It follows from Eq. (34) that SE oc n 2D (pf = 
\j2nn2D)- Electrons and holes may appear as charge 
carriers in graphene at T = when it is situated on 
a substrate and the electric field is applied (the elec- 
tric field effect). In this case, the concentration of the 
charge carriers is proportional to the applied voltage, 
n2D oc V g [1]. In the absence of a substrate, U2D = 
and 8E = and graphene in a vacuum remains a zero- 
gap semiconductor. 7 Thus, a substrate material with a 
sufficiently small n e ff value should be found. It could 
be, e.g., a substrate with regularly spaced pin holes 
or some metamaterials. The application of an electric 
field is unnecessary for the appearance of charge carri- 
ers in graphene in the case of the second-type contact 
of graphene with the substrate, i.e., when the e = 
level of graphene does not lie in the band gap of the 
substrate material. 

7. DISCUSSION 

The main qualitative result of this work is the pres- 
ence of the phase transition from a zero-gap semicon- 
ductor to a semimetal in the two-dimensional case at 
a* ~ 1. The zero-gap semiconductor phase is stable at 
a* < l, 8 which is an analog of the stability criterion for 
the system of three-dimensional non-relativistic fermions 



£tct(q,0) 

where £t D t (q,0) is the statistical limit of the total di- 
electric function. 

A semimetal-scmiconductor transition in disordered 
degenerate semiconductors was studied by Fradkin [14]. 
The HgTe (III-V semiconductors) and SnTe (IV-VI 
semiconductors) alloys, as well as the twodimensional 
graphite, i.e., graphene, were considered as the exam- 
ples of such systems. The transition appears due to 
the scattering of the charge carriers on a random po- 
tential and is caused by the carrier localization, but the 
Coulomb and spin-orbit interactions were disregarded 
in the model. In this work, the influence of the disorder 
in graphene on its transport properties is not consid- 
ered, but it is shown that the inclusion of the Coulomb 
interaction leads to a qualitatively similar result. 

7 There is also a purely technological difficulty in performing 
the experiment with graphene suspended in a vacuum: the pres- 
ence of a large number of defects such as vacancies and local 
ruptures is possible. It is not excluded that the graphene film 
is bent, creating a geometric potential. If all of these difficulties 
were overcome, electrons could be "sputtered" onto graphene, 
the necessity of the substrate as a source of the charge carri- 
ers disappears, and the semimetallic state of graphene becomes 
possible. 

8 When a* exceeds unity by an infinitesimal value, the system 
becomes unstable. However, at the exact equality a* = 1, E gs = 
should be expected and the system still remains stable. 
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Abrikosov and Bcneslavskii [15] considered the case 
of a Fermi point, a single-point touching of the con- 
duction and valence bands at Ef = 0. In this case, the 
Fermi energy does not enter the Green's function and 
the subsequent calculations were performed at pp = 0, 
which is essentially different from the problem consid- 
ered in this work, namely, graphene on a substrate 
with a non-zero concentration of charge carriers, i.e., 

8. CONCLUSIONS 

In this work, the ground-state energy of an elec- 
tron (hole) gas was calculated including the renormal- 
ization of the quantity u. The stability of a zero-gap 
semiconductor phase with respect to the transition to 
a semimetallic state was analyzed on the basis of the 
calculation. The condition of the transition appearance 
was obtained in the form a* ~ 1. Possible systems for 
the experimental observation were pointed out. 

I am deeply grateful to A. P. Silin for valuable advice 
during the course of this work and fruitful discussions 
of the results. 
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